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CO . The ionization history of the universe provides a major source of ambiguity in constraining cosmological 

, models using small angular scale microwave background anisotropics. To clarify these issues, we consider a 

, complete treatment of Compton scattering to second order, an approach which may be applicable to other 

Qv^ . astrophysical situations. We find that only the 0{v) Doppler effect, and the 0{vS) Vishniac effect are 

r~| . important for recent last scattering epochs; the C(f^) Doppler effect is not significant on any angular scale, 

^p^, and other higher-order effects are completely negligible. However the 0{v^) effect does lead to Compton- 

o 



y distortions, which although generally below current constraints, set an unavoidable minimum level in 
, reionization models. We consider the small-angle approximation for the Vishniac effect in detail, and show 

^ several improvements over previous treatments, particularly for low f^Q. For standard cold dark matter 



models, the effect of reionization is to redistribute the anisotropics to arcminute scales; late reionization 
leads to partially erased primary fluctuations and a secondary contribution of comparable magnitude. Using 
5^ , recent anisotropy limits from the ATCA experiment, we set new constraints on baryonic dark matter models. 

Stronger constraints are imposed (in second order) upon models with higher Hubble constant, steeper n, 
and higher density. These limits depend on the specific ionization history assumed, but the factor gained by 
lowering the ionization fraction is generally small, and may be tested by currently-planned experiments on 
arcminute scales. 

Mingled and merged, densely sprouting 
In the primeval mass, there is no shape 
Spreading and scattering, leaving no trail behind 
In the darkness of its depths there is no sound 
Chuang-tzu 
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I. INTRODUCTION 

Temperature fluctuations in the cosmic microwave background (CMB) are a direct probe of density 
perturbations at z ~ 1100. Therefore measurements or upper limits on such fluctuations can be used to 
constrain cosmological models for the evolution of structure. However, there is a possible loop-hole at small 
angular scales, since rcionization could make the last scattering epoch more recent, and erase the small-scale 
fluctuations, e.g., [1]. This possibility requires significant levels of ionization back to 2; ^ 10, to reach 
optical depth unity at z ^ 1100, which is certainly feasible in some models. With the recent detection 
of temperatur(^ anisotropics at large angular scales by Differential Microwave Radiometer (DMR) on the 
Cosmic Background Explorer (COBE) satellite [2], the constraining of cosmological models has entered a 
new phase of precision, with application to intermediate angular scale anisotropy experiments, e.g., [3]. It is 
becoming increasingly important to know how much leeway rcionization scenarios can give for small angular 
scale limits. In this spirit, we have attempted to systematically consider the effects of Compton scattering 
on CMB photons. 

Primary fluctuations arising from the recombination epoch are erased on scales smaller than that sub- 
tended by the horizon at the new last scattering epoch. However, secondary fluctuations will be generated 
on this new surface of last scattering, mainly due to Doppler shifts among the scatterers. These secondary 
fluctuations are 0(v), except that there is a partial cancellation of blucshifts and redshifts through a given 
overdense (or underdense) region. The secondary fluctuations were therefore thought to be small until Os- 
triker & Vishniac [4] pointed out that a second order contribution [of 0(vS), the so-called Vishniac term], 
which does not suffer from the cancellation, could be larger than the flrst order term. In practice the second 
order term is of the same order of magnitude as the degree-scale primary fluctuations which have been erased, 
but on arcminute scales and below. 

Use of reionization to avoid small-scale AT/T limits therefore depends on a calculation of these second- 
order contributions for the particular power spectrum and ionization history being considered. Detailed 
calculations for some models were performed by Vishniac [5] and Efstathiou [6]. However, the Vishniac 
term is only one of a number of possible second-order terms. Although it has been assumed that this term 
dominates over any others, this has not been demonstrated. Here we show that it is reasonable to neglect 
all other second order terms. Consideration of secondary anisotropics is therefore calculated to reasonable 
accuracy (in the linear regime) if the first order and Vishniac contributions to AT/T are evaluated. To arrive 
at this result, we set up in §11 a methodology which finds all second-order Compton scattering terms. Along 
the way, we uncover a number of physical effects which have not previously been described, e.g., minimal 
spectral distortions required in reionization scenarios independent of thermal history. Our general results, 
detailed in Appendices A and B, may also be of use for a wider range of problems. In §111, we discuss 
the second order contributions to anisotropy and significantly improve the approximations of Efstathiou [6] . 
Predictions for the cold dark matter (CDM) and baryonic dark matter with isocurvature fluctuations (BDM, 
also known as primordial baryon isocurvature) scenarios are computed in §IV. Recent limits on arcminute 
scale fluctuations from the Australian Telescope Compact Array (ATCA) [7] place strong constraints on 
BDM models. Other studies of the effects of reionization on the microwave background [5, 6] have tended 
to concentrate on the simplest case of a universe in which there has been no recombination. Here we have 
also considered the effects of more realistic ionization histories. 

II. THE BOLTZMANN EQUATION 
FOR SECOND ORDER COMPTON SCATTERING 

A. GENERAL FORMALISM 

The Boltzmann equation in general is given by 

dt ^ dt ^ dji dt ^ dpo dt ^ ' 

where / is the photon occupation number, 7; are the direction cosines for a photon of 4-momentum p, and 
the expression C{x,p) is the usual collision term. Latin indices range from 1-3, and we have employed 
the implicit summation convention. There are two distinct classes of second order Boltzmann equations 
that we might consider. The left hand side of equation (1) may be expanded to second order in metric 
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perturbations; whereas the right liand side, in the case of Compton scattering, may be expanded to second 
order in the energy transfer from coUisions. In Appendix A, we carry through a full derivation of all second 
order gravitational terms due to metric perturbations in the synchronous gauge. However, in this paper we 
are primarily interested in the effects of Compton scattering which are constrained by causality to manifest 
themselves on small scales where gravitational effects are unimportant. Thus we take the zeroth order 
approximation of equation (1) with respect to the metric fluctuations and find 

df df dx' 1 da df 

The third term on the left merely expresses the cosmological redshift of the photon energy, po oc a^^, where 
a{t) is the usual scale factor. In the homogeneous and isotropic limit, the second term in equation (2) 
vanishes. Spectral distortions in the; early universe are usually computed in this approximation (see part 
B). In §111, we calculate the effects of dropping the assumption of homogeneity. If, on the other hand, we 
are only concerned with temperature and not spectral distortions, we may integrate over momentum p. The 
third term may thus be eliminated since temperature and energy redshift in the same manner. 

The goal is to derive the collision term for Compton scattering, j{p) + e{q) j{p') + e(g'), to second 
order in the small energy transfer due to scattering. Note that we are performing calculations in the linear 
regime. Thus, the results are of interest only when the effects of the first order terms suffer cancellation, 
as in the case of a thick last scattering surface. The lowest order term is then of second order. Third and 
higher order effects are therefore negligible as long as the second order term is not cancelled. 

Our approach may be of more general interest since it provides a coherent framework for all Compton 
scattering effects, be they spectral distortions or anisotropies. In the proper limits, the equation derived 
below reduces to familiar equations and effects, e.g., the Kompaneets equation (Sunyaev-Zel'dovich effect), 
linear Doppler equation (Vishniac effect). Furthermore, new truly second order effects such as the O(w^) 
quadratic Doppler effect are obtained. These effects predict distortions well below the observational limits 
today. In principle, however, the fact that they may mix both anisotropies and spectral distortions makes 
them distinguishable. 

We make the following assumptions in deriving the equations: (1) the Thomson limit applies, i.e., 
the fractional energy transfer 6p/p <i; 1 in the rest frame of the background radiation; (2) the radiation 
is unpolarized and remains so; (3) the density of electrons is low so that Pauli suppression terms may be 
ignored; and (4) the electron distribution is thermal about some bulk flow velocity v. Approximations (1), 
(3), and (4) are valid for most situations of cosmological interest. The approximation regarding polarization 
could be dropped to give coupled equations for perturbations in the total and polarized components [8]. 
Polarization perturbations are typically an order of magnitude smaller than temperature distortions [9]. 
Thus the contribution from polarization to the evolution of temperature distortions is a small effect, as 
the microwave background never generates a significant polarization. For calculations on the generation of 
polarization in second order theory, see [6]. 

The collision term may in general be expressed as [10] 

^^""'P^ = 2^ / DqDq'Dp'{2nr6^^\p + q-p' -q')\M\' 

X {5(x, q')/(x, pO [1 + /(x, p)] - 5(x, q)/(x, p) [1 + /(x, p')]} , 

where |Mp is the Lorentz invariant matrix element, /(x, p) is the photon distribution function, ^(x, q) is 
the electron distribution function and 

d^q 

" (2^)32i?(q) 

is the Lorentz invariant phase space element. The terms in (3) which contain the distribution functions are 
just the contributions from scattering into and out of the momentum state p including stimulated emission 
effects. 
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We will assume that the electrons are thermally distributed about some bulk flow velocity v, 

i?(x,q) = (27r)=^ne(27rmTe)-3/2exp| ^^LZ^^I^^l , (4) 

where m is the electron mass, and we employ units with c = h = k = 1. Expressed in the rest frame of the 
electron, the matrix element for Compton scattering summed over polarization is given by [11] 

(5) 

where the tilde denotes quantities in the rest frame of the electron, a is the fine structure constant, and 
cos/3 = 7 • 7' is the scattering angle. Note that here and below wc define p = po = |p|. Of course, the 
matrix element must be expressed in terms of the corresponding quantities in the frame of the radiation for 
calculation purposes [see equation (B-1)]. 

The result of integrating over the electron momenta can be written 

C^P'— j [^0 + Cp/m + Cy + Cyy + CT^/m + Cyp/m + C(^p/my] , (6) 

where we have kept terms to second order in Sp/p. The explicit expressions for the quantities in equation (6) 
are given in equations (B 4) and (B 5) of Appendix B and arc discussed in turn below. This equation may 
be considered as the source equation for all first and second order Compton scattering effects. However, in 
most cases of interest only a few of these terms will ever contribute. We have therefore written the collision 
equation implicitly and expressed it in terms of the sealing behavior of the contributing elements. 

The expansion of the Compton collision term to second order in 6p/p has actually involved several 
small quantities. It is worthwhile to compare these terms. The quantity Te/m characterizes the kinetic 
energy of the electrons and is to be compared with p/m or essentially T/m where T is the temperature 
of the photons. Before a redshift 8.0(r2o/i^)^^^a;^^^^, where Xe is the ionization fraction (this corresponds 
to 2; > 500(f2B/i^)^/^ for standard recombination), the tight coupling between photons and electrons via 
Compton scattering requires these two temperatures to be comparable. At lower redshifts, it is possible 
that Te T, which produces distortions in the radiation via the Sunyaev-Zel'dovich (SZ) effect as discussed 
further below. For the reionization scenarios which we consider, these two temperatures will typically still 
be comparable at last scattering, ~ 50 with T/m ~ 5 x 10~"'^''(1 + z*). The SZ effect, however, will play 
a role at lower redshifts even in these scenarios. Note that the term Tg/m may also be thought of as the 
average thermal velocity squared (wj^) = 3Te/m. This is to be compared with the bulk velocity squared 
and will depend on the specific means of ionization. The bulk velocity is related to the fractional overdensity 
S = 6p/ p hy the continuity equation. On scales much smaller than the horizon, v <^ 6 in linear theory. 

It is appropriate at this point to examine the physical significance and qualitative features of each of 
the sources in the collision term of equation (6): 

1. Co- Anisotropy Suppression. 

In the absence of electron motion, there is no preferred direction. Thus to zeroth order, scattering 
makes the radiation distribution more isotropic. The Co term given by equation (B-4) of Appendix B 
indeed equalizes the distribution function over all directions via scattering into and out of a given mode. It 
is therefore responsible for the suppression of primordial anisotropics by reionization. It can also significantly 
affect the regeneration of inhomogeneities when ordinary contributions are suppressed, as in the case of thick 
last scattering surfaces (see §IIIC). Recall that an inhomogeneity on the last scattering surface becomes an 
anisotropy in the CMB today. 



|M|2 = 2(47r)2a2 
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2. Cy and C^^: Linear and Quadratic Doppler Effect. 

Aside from the small electron recoil (see 3.), the kinematics of Thomson scattering require that no 
energy be transferred in the rest frame of the electron i.e., p' = p. Nevertheless, the transformation into the 
background frame induces a Doppler shift, so that 

7 = 1 ^ J.'y - 1 = V • (V - 7) + (V • y)v • (V - 7) + 0{v^). (7) 

Notice that in addition to the usual linear term in v, there is also a term quadratic in v. 

To gain physical insight into those Doppler ofFccts, we c;an approximate; the photons as isotropic and 
neglect the angular dependence of Thomson scattering, (we will postpone discussion of the precise effects 
until §111). Averaging over the incoming direction 7, we obtain 

^|^:.V7' + (v-7')'. (8) 

Thus the linear Doppler effect introduces an energy shift whose sign depends on direction; whereas the 
quadratic Doppler effect, which is positive definite, always gives rise to a blueshift of the photons. 

Now let us consider the case where there arc many scattering blobs so that the directions of the electron 
velocities are in effect randomized. In this case, the net linear effect vanishes and 




i.e., since rcdshifts and blueshifts cancel, there is no net transfer of energy to the photons to first order in 
V. The quadratic Doppler term however represents a net energy gain of 

A = f = ^(A (10) 

since the energy density £ oc T^, and the Doppler shift of a blackbody is a blackbody with T = To{l + 5p/p). 
Here, we assume for simplicity that all the photons in the spectrum scattered once. 

These effects are not wholly equivalent to a uniform Doppler shift. In averaging over angles above, we 
have really superimposed many Doppler shifts for individual scattering events. Therefore the resulting spec- 
trum is a superposition of blackbodies with a range of temperatures AT/T = 0{v). Zel'dovich, Illarionov, 
& Sunyaev [12] have shown that this sort of superposition leads to spectral distortions of the Compton-y 
type with y = O(v^). Thus, spectral distortions have a quadratic dependence on v and may be considered 
as part of the quadratic Doppler effect. 

In summary, the linear Doppler effect is primarily cancelled after many scatterings. Nevertheless, some 
residual effects remain due to the evolution of the electron velocities and densities during the last scattering 
epoch [13] and the increased probability of scattering in an overdense region [4] . The quadratic Doppler effect 
leads to an average net energy increase of A 2± 4/3t(w^), and a spectral distortion of the Compton-y type. 
Here r = / neaxdt is the optical depth, which gives the fraction of photons scattered if t 4C 1. We will now 
see that the quadratic Doppler effect for bulk flows is in fact entirely equivalent to the Sunyaev-Zel'dovich 
effect for thermal motions. 
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3. Ct^Itu and Cp/m' Thermal Doppler Effect and Recoil. 

Of course, we have artificially separated out the bulk and thermal components of the electron velocity. 
The thermal velocity leads to a quadratic Doppler effect exactly as described above if we make the replace- 
ment {v"^) (vth) = ^Te/m [e.g., in equation (10)]. For an isotropic distribution of photons, this leads to the 
familiar Sunyaev-Zel'dovich (SZ) effect [14]. The SZ effect can therefore be understood as the second order 
spectral distortion and energy transfer due to the superposition of Doppler shifts from individual scattering 
events off electrons in thermal motion. This can also be naturally interpreted macrophysically: hot electrons 
transfer energy to the photons via Compton scattering. Spectral distortions result since low energy photons 
are shifted upward in frequency, leading to the Rayleigh- Jeans depletion and the Wien tail enhancement 
characteristic of Compton-?/ distortions. The fractional energy change due to scattering is 

A % ^ ^tKJ. (11) 

Note that the thermal and quadratic Doppler effects are wholly equivalent in this sense. 

Even though this is a net energy gain, the effective temperature in the Rayleigh-Jeans region is sup- 
pressed due to photon depletion (up-scattering), 

Ar\ A 

— =-2yc.--. (12) 

This is to be distinguished from the linear Doppler effect which provides a frequency independent shift in 
temperature given by AT/T ~ A/4 (c./. §§IIA2). 

If the photons have energies comparable to the electrons [i.e., the electron and photon temperatures are 
nearly equal), there is also a significant correction due to the recoil of the electron. The scattering kinematics 
tell us that 



l-f-—(l- cos/3) 
m 



(13) 



Thus to lowest order, the recoil effects are 0{p/m). Together with the thermal Doppler effect, these terms 
form the familiar Kompaneets equation [15] in the limit where the radiation is isotropic (as described in part 
B). The combination of thermal Doppler shift and recoil drive the photons to a Bose-Einstein distribution of 
temperature Tg. A blackbody distribution cannot generally be established since Compton scattering requires 
conservation of the photon number. 

4. C^p/to and C(p/„)2: Higher Order Recoil Effects. 

These terms represent the next order in corrections due to the recoil effect. In almost all cases, they are 
entirely negligible. Specifically, for the CMB, when electron bulk flow effects become significant (p/m) < 
0{v'^). Furthermore, since there is no cancellation in the Cp/^ term, (7(^/^)2 will never produce the dominant 
effect. We will hereafter drop these terms in our consideration. 



B. THE ISOTROPIC LIMIT AND SPECTRAL DISTORTIONS 

The most interesting case to consider is when the radiation is nearly isotropic. For an initially isotropic 
radiation field, only the Doppler effects, which have a preferred direction along the velocity of the electron, 
can generate an anisotropy. However, as we have seen, the first order contributions due to the linear Doppler 
effect tend to cancel. Furthermore, to lowest order, scattering reduces anisotropics. It is therefore reasonable 
to consider the case where the anisotropy itself is at most second order. We will justify this more fully in 
§111. The general results of Appendix B can be employed to study the effects of Compton scattering on 
highly anisotropic radiation. For example, equation (B-4) contains the anisotropic generalization of the 
Kompaneets equation which may be useful in other astrophysical settings. 
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Under the additional assumption of near isotropy, the colhsion term becomes more tractable and trans- 
parent: 



C(x,p) = neaT< 



/0-/+^P2(m)/2 



20 20^' ^ 



mp^ dp 



P^Te^^ + /(! + /) 



(14) 



where 



fi 



(15) 



are the Legendre moments of the distribution function with respect to the axis 0:3, = 7 ■ X3. We assume 
that for plane wave situations there is an azimuthal symmetry about the direction of the wavevector X3 since 
to lowest order the source term possesses this symmetry (see §IV). 

The first term in equation (14) just represents the zeroth order suppression and can be found in the 
standard treatments in the literature, e.g., [16]. In the limit of complete isotropy of the radiation field, 
it vanishes identically. The second and third terms represent the linear and quadratic Doppler effects 
respectively. The final term is the usual Kompaneets equation. 

Notice that in the limit of many scattering blobs, we can average over the direction of the electron 
velocity. The first order linear Doppler effect primarily cancels in this case. Assuming complete isotropy and 
ignoring the residual contribution of the linear Doppler effect, we can reduce equation (14) to 



C(x,p) = neCrT< 



1 d 




1 d 




1 Z p^ dp 


dp 


mp^ dp 





T^^f +/(! + /) 



(16) 



Under the replacement (w^^) = "iTe^jm w^, the SZ (thermal Doppler) portion of the Kompaneets equation 
and quadratic Doppler equation have the same form. We also regain the factor of 1/3 described in (9). Thus, 
spectral distortions due to bulk flow have exactly the same form as SZ distortions and can be characterized 
by the Compton-y parameter given in its full form by 



Hear 



T 



dt. 



(17) 



A full analysis of the Boltzmann equation analogous to that performed in §111 yields this result as well. 
The appearance of the photon temperature T in equation (17) is due to the recoil terms in the Kompaneets 
equation [c.f. equation (11)]. 

As was noted by Zel'dovich & Sunyaev [17] the observational constraints on the Compton-y parameter in 
conjunction with other spectral limits can rule out some reionization histories independently of the detailed 
thermal history of the models. However, Bartlett and Stebbins [18] show that for a low ils universe as 
implied by nucleosynthesis, no such model independent constraints exist. Furthermore, for the case where 
the electron temperature and photon temperature are equal, there is no SZ effect as one can clearly see from 
equation (17). Therefore, the bulk flow effect places a lower limit on spectral distortions due to any given 
reionization history. 

For example in an Qq = ^ universe with full ionization up to a redshift Zi and Tg = T, 



y ~ 4.1 X 10-2(Ob/i) [(1 + ^i)V2 _ 1] (^;2^, 



(18) 



where (vq) is the (spatially) averaged square of the velocity today. Here and throughout the paper Qb refers 
to the fraction of the critical density of the intergalactic medium in baryons; it does not include baryonic 
matter in compact objects. We have assumed linear growth of velocities in an fio = 1 universe, 



v{k,r]) =i^ = ^ OCT] [^^0 = 1], 



(19) 
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where the dot clenotcs a derivative with respect to coiifornial time rj = J dt/a. Given a power spectrum P{k) 
for the baryon density fluctuations, we can compute the average velocity as 




assuming ergodicity. Here Vx is the volume in which the universe is assumed to be periodic. The result 
for the standard cold dark matter (CDM) scenario, normalized to the COBE measurement [2], is therefore 
vo,rms = 1100 km/s assuming that the baryons follow the dark matter distribution (see §V for details). Note 
that this value is rather large, since COBE requires a bias factor 6 1 which produces excessive small scale 
velocities. For comparison purposes, the sun's velocity with respect to the CMB is 365 ± 18 km/s whereas 
the velocity of the local group is 622 ± 20 km/ s [19] . The upper limit on the Compton-y distortion from 
the Far Infrared Absolute Spectrophotometer (FIRAS) [20] is y < 2.5 x 10~^. Therefore it is appropriate to 
reexpress equation (18) as 



where z„,,„.t is the maximum redshift that the electrons can be fully ionized. 

However our assumption that the electron density follows the dark matter is invalid at high redshifts. At 
redshifts z > Zdrag — 160(r^o/i^)^/^a;^^^^, Compton drag on the electrons and protons prevent matter from 
falling into the dark matter wells. Thus, (21) will not give a constraint if z,,,,,,,. > Zdrag- Correspondingly, 
the prediction for CDM is y{zdrag) — 3 x 10~^, almost two orders of magnitude below the current limits. 
Almost certainly foreground contamination from hot clusters will make this effect unobservable. Estimates 
employing cluster modeling yield y ~ 2 — 6 x 10~^, but are extremely sensitive to the normalization of the 
power spectrum and cluster dynamics [21]. Therefore in flat fio = 1 models, the quadratic Doppler effect 
does not yield a measurable average (isotropic) Compton-y parameter. In §IIIE, we calculate the fluctuations 
in the Compton-y parameter at arcminutc scales and show that it too is negligible. Note that the isotropic 
Compton-y distortion is an integral over optical depth dr, i.e., it has a contribution from all scatterings. 
On the other hand, correlations {i.e., temperature fluctuations) come from an integral over the visibility 
function exp(— T)(iT, i.e., they are dominated by the last scattering event. 

The v"^ Compton-y distortion is non-negligible in the case of open universes. In an open universe, 
velocities grow more slowly so that they are still reasonably large at the scattering epoch. Furthermore last 
scattering happens at lower redshifts for baryon dominated open models. In models with relatively low CIb, 
predictions for the quadratic Doppler effect are comparable to estimates of the SZ effect, although still below 
the present constraints; higher Qb models with an early photoionized IGM can already be ruled out from 
their larger thermal y-distortion (see e.g., [22]). 

It is worthwhile to note that such quadratic Doppler distortions are the minimum required in the case 
of reionization and unlike the SZ eff'ect do not suffer from uncertainties in cluster dynamics and ionization 
mechanisms. 

III. THE BRIGHTNESS EQUATION AND SMALL SCALE ANISOTROPIES 

A. GENERAL ARGUMENTS AND DEFINITIONS 

In this section, we will extend, refine and unify techniques developed in refs. [5, 6, 13]. We employ the 
true angular dependence of Compton scattering and treat the 0{v'^) quadratic Doppler effect as well as the 
0{v) and 0{v6) effects. Furthermore, we find that for Qq < 1 models, there are significant corrections to 
the formalism of Efstathiou [6] . 

Having shown that spectral distortions are minimal and simply described by the Compton-y parameter, 
we will concentrate on temperature anisotropics, i.e., we shall be no longer be concerned with frequency 




(21) 
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dependence. Therefore let ns integrate equation (14) over the spectriim to obtain the evohition equation for 
the fractional energy perturbation to the spectrum A = 5£/£. This yields the so called brightness equation 



A + ^idi/^ = nearall + <5(x)]{Ao - A + ^P2(/u)A2 
+ 47 • V - + 7(7 • v)2|, 



(22) 



where the dot denotes differentiation with respect to conformal time, and we have assumed that the cluster 
SZ effect is unimportant at the last scattering epoch. The average SZ effect will produce an isotropic 
Compton-y distortion and have no effect on small scale anisotropics. 

Here, we separate the average electron density from the spatially dependent perturbation as 



rie = a;e(?7)(l - yp/2)ns(x, t?) = ne{r])[l + S{x, r?)]. 



(23) 



where Yp is the primordial mass fraction in helium, Xe is the ionization fraction (assumed to be independent 
of position), is the baryon number density and S is the fractional overdensity in baryons. The spatial 
dependence of the electron density distribution couples with the ordinary source terms to provide an addi- 
tional effective source. The coupling of density perturbations to the velocity field is in fact responsible for 

the Vishniac effect. 

Parenthetically, note that our definition of A does not exactly conform with the standard convention [6] , 
where A = 5f{T/A){df/dT) ~ ^{T^s{v) - T)/T, with / a Planck distribution at the average temperature 
T and Toff (p) the effective temperature of the spectrum as a function of frequency. Under this definition, A 
retains spectral information. We define A as the perturbation in the total energy density, which is explicitly 
independent of frequency. The two definitions coincide for a uniform shift in the temperature {e.g., the linear 
Doppler effect), which is the situation usually considered. Care should therefore be taken in associating our 
A with temperature distortions (c./. equation (12)). 

We can rewrite equation (22) implicitly as 



A + 7,a,A + f(r?)A = f(??)5(x, 7, r?). 



(24) 



where the conformal time derivative of the optical depth, f(?7) = neCxa, can be interpreted as the probability 
of scattering in the interval drj. Note that the source term 5(x, 7,77) includes the anisotropy suppression 
(Ao) part. The other source terms, as we shall see, generate a small contribution to the isotropic fluctuations 
by scattering out of a given mode, thereby partially escaping the cancellation effect. The isotropic part then 
contributes by scattering back into the line of sight. This of course does not increase the anisotropy, only 
the inhomogeneity of the photon distribution. 

To solve this equation, we transform to its Fourier space analogue, 



A(k,7,7?) = -|^ j d^a;A(x,7,r/)exp(-ik-x) 
which has the well known solution [6], 



A(k,7,?7)= / S'(k,7,r7')5(»?.'7')exp[zA:/z(r7' - ?7)]c^??', 

'0 



where /U = 7 • k/fc and 



3(77, 77') = f(?7)exp 



dr]"f{r]") 



rj' 



— exp[-T(77,77')] 



(25) 



(26) 



(27) 



is the so-called visibility function [note that g{r], rf) = f{r])]. For the case of full ionization where Xe = 1, 
t{z) = ^^(1 - ^p/2)Sf [2 - 3^^o + (1 + ^oz f^^oz + 3^0 - 2)], 



(28) 
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(where rup is the mass of the proton) or alternatively 



where TZ = Hq is the curvature scale today which merely translates our convention for conformal 

time to the usual development angle ^ = ri/TZ. For Qq = 1, the expression in braces should be replaced by 
iiVo/v? - 1}. 

We are primarily interested in the solution at small wavelengths, i.e., k 6r] ^ 1 where dr] is the conformal 
time "thickness" of the last scattering surface, located at 77*. The precise definition of 77* is somewhat 
arbitrary, e.g., 77* could be defined as the epoch at which optical depth becomes unity (we will employ this 
definition unless otherwise stated). For late last scattering St] c±r]^; it is in this sense that we mean the last 
scattering surface for reionization scenarios is "thick." 

Neglecting the evolution of the source term and the optical depth, the integral over the plane wave 
for /X 7^ in equation (26) will in general cancel except for a small portion of the order of a wavelength. 
Therefore the resulting brightness fluctuations will be of order 

A(k,7,r7)~5(k,7,77)rA(r7), ^x^0, (30) 

where t\ = {2TT/k)f{r]) is the optical depth across one wavelength. If the perturbation has many wavelengths 
over the last scattering surface, ta <C 1 since the optical depth across the whole surface is of order unity. 
Therefore, the brightness fluctuations arc suppressed as compared with the source terms. 

The n = mode escapes this cancellation. Thus the brightness fluctuation has a significant value only 
if the wavevector k is perpendicular to the direction of the photon momentum, 7. The interpretation of this 
statement is simply that if the wavevector were not strictly perpendicular to the photon momentum vec;tor, 
photons in a given direction would have scattered off both crests and troughs of the perturbation, yielding 
net cancellation. 

In the case of the first order (linear Doppler) term, S'(k, 7, ry) = 47- v = 4/iw which vanishes for the /i = 
mode. Therefore, as pointed out by Kaiser [13], the first order term is primarily cancelled up to a factor which 
involves the time evolution of the perturbation. It is important however to note that the general cancellation 
argument holds for any source term, while near perfect cancellation will occur for S'(k, /i = 0, 77) = 0. 

Having determined which perturbation modes give large contributions to brightness fluctuations, we 
need to translate this information into the angular correlation function of temperature fluctuations on the 
sky today. This is denoted by 

/AT AT \ 
C{e,a) = (^ — {^,) — M)^ , (31) 

where the average is performed with fixed beam throw — cos^^(7i • 72) and with Gaussian beam-width a 
[C{6, a) is not to be confused with the Compton scattering collision terms C(x, p)]. First we must relate the 
brightness distortion A to a fractional temperature perturbation. For linear Doppler effects which change the 
temperature uniformly across the spectrum and leave the spectrum a blackbody, A ~ AAT/T. Quadratic 
and thermal Doppler effects create Compton- 7/ distortions such that A ~ —2{AT/T)rj [see equation (12)]. 
To keep the expression applicable to both sorts of distortions, we will write 

^ = ^A, (32) 

where K = lior the linear Doppler effect and K = —2 for quadratic and thermal Doppler terms (Compton-y 
distortions). 

We will use the small scale approximation given by Doroshkevich, Zel'dovich, & Sunyaev [23], 
^^^'"^ Jo ^'^^^7 ^rf/x|A(fc,M,r?)|'Jo [kRr,e{l - exp [-(fci?^a)2(l - m')] , (33) 
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where i?,, translates comoving distance at epoch r] to angle on the sky. It is given by, 



or equivalently 

- I 7esinh[(,yo - V)/T^l ^0 < 1- ^ ' 

Efstathiou [6] takes the asymptotic form of equation (35) for rjo :$> rj [or (l^o-^)^^^ ^ !]• This is not a good 
approximation for Oq < 1 models in which the surface of last scattering was relatively recent. For example, 
if no = 0.2 and h = 0.8, then ~ 20. 

For Compton-y distortions, the correlation function [eqn. (33)] will in principle measure the uniform 
background ^/-distortion, estimated by equation (18), as well as fluctuations in the background. If no spectral 
information is obtained, it is not possible to observe the effects of the uniform distortion. However, in practice, 
this is unlikely to be a concern; in small-scale experiments the uniform background is always subtracted out. 
For example, many experiments employ beam switching to minimize noise, e.g., a double beam experiment 
measures C(0,(t) — C{9,<t). Whereas C{0,a) is sensitive to fluctuations from all scales greater than cr, 
C{0,a) — C{0,a) is effectively only sensitive to angular scales between a and 9. Even maps generated by 
interferometry are only sensitive to fluctuations below the size of the primary beam Opb- 

The cancellation arguments, which say that |A(fc, fi, r])\ is sharply peaked at /U = 0, allow us to set /U = 
in the slowly varying terms, so that 

POO 

C{e,a)= Q{k)Jo{kRr^e)exp[-{kRr,af]dk/k, (36) 
Jo 

where 

Q{k) = £^ \A{k,t,,v)\'dn (37) 

expresses the power per logarithmic interval of the fluctuations. 

Note that technically speaking since the square of the brightness distortion enters into (37), we need 
to worry about cross terms between first and third order processes. However, they arc suppressed for two 
reasons. Firstly, we arc only interested in second order effects when the first order term is suppressed. The 
cross terms are thereby also suppressed. Secondly, as we shall see in part B, for the linear Dopplcr effect, 
the third order terms that couple with the first order term are also themselves suppressed. We will denote 
the order of the term by superscripts. For example, 

A = A^^) A(2) . . . , 

3 = 3W +S(^) + ..., (38) 
v = v«+v(2) + ... 



where superscript (1) refers to the values given by linear theory and (2) refers to the mildly nonlinear next 
order correction, etc. 

Finally, it should bo mentioned that equation (33) is only strictly valid for a thin last scattering surface. 
It assumes that there is a one to one correspondence between angles and length scales at last scattering 
given by equation (35). It also assumes that the phases of neighboring modes have been correlated by the 
scattering. We expect errors from these approximation to be of order 6r]/{r]o — r?*) — (~ 0.2 for 

realistic reionization scenarios). 
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B. GEOMETRICAL CONSIDERATIONS 



The physical picture becomes much clearer if we examine the geometrical aspects of the linear Doppler 
effects. Evolution under the linear Doppler and isotropic source terms can be expressed as 

A + 7i5iA + f(77)A = f(?7)[Ao + ^A2P2(m) + 47 • j(x)], (39) 

where 

j(x) = [1 + (5(x)]v(x) = v(i) + (5(i)v(i) + v(2) + . . . (40) 

is the so-called normalized matter current introduced by Vishniac [5]. As discussed in part A and shown 

explicitly below in part C, the source terms involving A itself arc small compared to the Doppler terms, and 
thus higher order terms in A have been dropped. We can reinterpret equation (39) as representing scattering 
off a uniformly dense medium of electrons with velocity j(x). The increase in probability of scattering off 
overdense regions is then accounted for by a rcscaling of the velocity field. 

For the linear Doppler terms, we can also think of the cancellation argument [see discussion following 
equation (30)] as follows. If the matter current has zero vorticity, the contribution to the brightness fluctu- 
ation is approximately the line integral of the gradient of a field. The contribution is thus cancelled apart 
from contributions from the end points. This is equivalent to the cancellation argument: if V x j(x) = 
then j(k) || k and thus S{k, /x, ry) oc /x ~ 0. 

Now let us consider the matter current vorticity: Vxj(x) = [1 + (5(x)] [V x v(x)] — v(x) x V(5(x). Since 
gravitational perturbations of a pressureless fluid are irrotational, V x v(x) = to all orders. However 
this is not true of the second term. The vorticity of the (vS) term (and in second order theory this term 
alone) does not vanish. Although the velocity field is the gradient of a potential, it is not the gradient of 
the local density field. Small scale overdensities can move coherently with some peculiar motion along the 
line of sight. In these regions, there is an increased probability of scattering and inducing a Doppler shift 
in the background radiation temperature. This will imprint an anisotropy on the CMB at the scale of the 
scattering blobs. Note that this is an inherently non-linear process since we are coupling a large scale (the 
scale of the velocity field) with a small scale (the scale of the overdensity). Recall that the Fourier transform 
of a product is a convolution, and 

poo 

j„5(k)=v*,5= / d^k'v{k')6{k-k'), 
Jo 

which shows the mode coupling explicitly. Since v(k') || k', this gives a net non-zero contribution perpendic- 
ular to k. In other words, we can have plane waves of the vector field j(x), with wave vector perpendicular 
to the line of sight, but with components parallel to the line of sight. This is exactly as required above: there 
is no cancellation between crests and troughs for the mode k _L 7, but the Doppler effect does not vanish 
since 7 • j 7^ 0. 

It is also useful to note that the perpendicular and parallel components add in quadrature, i.e., the cross 
terms between the Vishniac term and other second order terms cancel when we calculate |A(A:, n, ry)^, as we 
now show. In all cases of interest, the source terms can be separated into a term dependent on conformal 
time and terms dependent on the mode coupling, 

rj.giv, v')S{k, 7, rj') = A{krj.)G{r = 7?7ry., r?) ^ F(7, k, k') S^'^ (k', rj*) S^'^ (k - k', ,?.). (41) 

k' 

where we have written the conformal time ratio r?'/^* ^ Therefore employing equation (26), we obtain 
the solution 

A(2) (k, 7, rj) A{kr].)G{kfiT],)eM-iknv) Yl ^(^' ^^^^ C^'' ^^^^ 0^ ' ^*)' (^2) 

k' 
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where we have approximated the integral over confornial time as a Fourier transform by taking r/ — > oo in 
both the Umits of the integral and in G, i.e., 



I G{r,r])exp{ikfj,r]t,r)dr / G(r, 77 = oo)exp(ifc/ur/*r)rfr 

Jo Jo (43) 

Hence G{k^iri^) is the Fourier transform of G{r,ri = 00) with r = rf'/rj^ as the Fomier conjugate of kfir]^. 
Note that the quadratic Doppler term also obeys equation (42) as we shall see in part E. On the other hand, 
the linear Doppler terms have the additional property that F(7, k, k') — 7 • D(k, k') with D(k, k') only a 
vector function of k and k'. Choosing our coordinate system so that k || X3 and 7 lies in the Xi — X3 plane, 
we can rewrite 

F(7, k, k') = a\^'> {k, k', cose)j ■ k + a*,^^ (fc, k', cos6')7 • k' 

' (44) 
= a\^' {k,k' ,cose)kn + a^l' {k,k' ,cose){l - ii^y/'^k'sm0cos(l), 

where k^ = fc'(sin6'cos(/>, sin^sin(/), 0) is the component of k' perpendicular to k, and a^^\ are arbitrary 

functions independent of the azimuthal angle (j). Thus we can separate the contributions into A^^^ = A||^^ + 

A^^ which are due to j || k and j -L k. 
Cross terms such as 

A||(')Af ocF|f(7,k,k')Fi(7,k,k') \5^'\k' , r],)f \S^'\k--k' ,r],)f k'^dk' smO d9dct> 

r2n (45) 



oc / cos 
Jo 



then vanish under integration over the azimuthal angle (p. Here we have used the random phase hypothesis 
to eliminate one of the sums over modes. Of course, cross-terms between first and second order vanish under 
this assumption as well. 

Now consider the mixed first and third order term. This can be expressed as 



A*(i)a(3) oc j \a\p{k, k', cos6')fc/i + af\k, k', cos6')(l - n'^y/'^k'smecoscj)^ dep. (46) 

The perpendicular part vanishes under the integration over azimuthal angle. However, the parallel part (for 
the linear Doppler effect) is suppressed, as we can see, by the additional factor of fi in equation (46). Hence 
mixed first and third order terms arc not only suppressed by the smallness of the first order term but also 
by a corresponding suppression in the third order portion. These contributions are thus entirely negligible. 

Since all relevant contributions add in quadrature, we can calculate them independently. In part C, we 
will calculate the first order Doppler and isotropic effects; in part D, the Vishniac effect; and in part E the 
quadratic Doppler effect. 



13 



C. FIRST ORDER DOPPLER AND ISOTROPIC EFFECTS 

In Fourier transform space, the solution to the residual first order effect is given by 



AW(k,M,,?) 



giv, v')exp[ikij{r]' - r])]dr]' . (47) 



Here we have retained the A contributions to the source in the first order integral despite the suppression 
discussed in section part A, since there is a corresponding suppression of the first order linear Doppler source 
term that makes these two comparable. Note furthermore that for the first order effect the source term is 
only a function of 7 • k = fc/x and not of 7 in general. 

Taking the zeroth moment with respect to fi of equation (47), we obtain 



A(^)(k,7?) = J\rj'g{v,v') [^i'^MHv - v')] - 4iv^'^ji[k{v - v')] + l^i'^MKv - V)]} , 



(48) 



where we have employed the identity 

1 

3n{z) = -{-iT j ^ex^p{ifj,z)Pn{lJ-)dn. (49) 

Due to the oscillation of the spherical bessel functions, in each case the integral only has strong contributions 
around k{'q — r]') < 1. Since kr]* ^ 1, this implies ~ ry' (and of course, this is just the cancellation argument 

in a different guise). Therefore we can take the slowly varying quantities g, A^q\ and v^^'^ out of the integral. 
As we have discussed in part A, the first term and third term are suppressed by a factor of order the optical 
depth across one wavelength of the perturbation. We can perform the integral over the second term in 
equation (48) by noting 

J, 2 r[i(n + 2)]' ^^^^ 

and therefore 

A«(k,,).-l!^^f(,). (51) 



Taking the second moment of equation (47), ignoring the A terms, and extracting the slowly varying 
quantities, we obtain 

A^^\k,rj):^Aiv^'\Kn)f{v)J\v'[ps[k{v-v')] - ^ii [M^ -??)]} • (52) 

However, employing equation (50), we see that A2^'' vanishes. The quadrupole term can thus be neglected 
as compared with the other sources. 

Now let us separate the time dependence of the perturbations. If we express the growth of density 
perturbations in linear theory by 

5W(k,r?)=I?(77,r7*)5W(k,7?*), (53) 



then the growth of velocity perturbations is 

' S^'\KV*), (54) 



^;W(k,r?)=ip 



where Vir],!]^,) = D{ri)/D{r]^) is the ratio of the growth factors [24]: 

'^^'^^°^{3sinh(??/7^)[sinh(7^/7^)-(??/7^)]/[cosh(77/7^)-l]'-2, nl<i. ^^^^ 
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Hence, for fio = 1 or rj/Tl <C 1 (alternatively (fio^:)"'^''^ ^ 1), 'D{rj,rj^,) = {rj/'q^)'^. Efstathiou [6] employs 
this approximation for ^Iq < 1 models. Again this is not a good approximation due to the fact that last 
scattering is so recent in these models. We use the correct growth factors here. Equation (47) can then be 
written as 



A«(k,/x,r?) 



8r5> -^'(k. '/*) 



''"'i' f {^^^''''''*^} [^(V) + iM5(r/,V)e'"^'''d7y'. (56) 



Let us check that this solution is consistent with the assumption of small Aj^^ . Taking the second moment 
of equation (56) and employing (49), we obtain 



4'^(k,7?) 



)(i)(k,r;.) TT d 
{kr]^y 4:k dr] 



(57) 



We see that A2^''/Aq^'' ~ 0(1/^77*) <IC 1 and Aj^^ is indeed small compared with the other source terms. 
The integral over the conformal time in equation (56) is again approximately a Fourier transform [13] 

of 



1 



dr 



d_ 
dr 



dr 



'D{rr]*,rj*) 



where the transform is defined in equation (43). Hence 



|A«(fc,M,ry)|^ = ^^glJi^|G.(fcM^.)^ 



(58) 



(59) 



where we have implicitly used the ergodic and random phase hypotheses. The power per logarithmic interval 
Q{k) for the first order Doppler effect is therefore 



(60) 



where 



v,v = 7r / Gy{kiJ.r]*)d{kiJ,r]*) 

J-kn, 



(61) 



\Gv{r,r]o)\'^dr, 



and we have used Parseval's theorem to approximate For full ionization with fio = 1 and r?* -C r?o, 
Iri.v — 6.38. Efstathiou [6] employs this value in cases where the approximation is not very good. We shall 
see in part E that this can lead to a significant correction. 

One might worry about the appearance of rj* in the expression for Q{k), since there is a certain arbitrari- 
ness in the definition of the epoch of last scattering, e.g., Efstathiou [6] picks t(t/*) = 1, whereas we might 
equally well have chosen the peak of the visibility function. Nevertheless, inspection of the growth factors 
in P{k,r]i,) and shows that all dependence on r]^ vanishes, as it must, since it has only been introduced 
as a convenient normalization epoch. 
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D. THE VISHNIAC EFFECT 

Contributions to the brightness fluctuation due to the Vishniac term (x, 7,7/) 

may be expressed as 



= 45W(x)[7.v(x)], 



where 



Gvs{r,r]) = 



^;5(7,k,k') = 7' 



dr 



'D{rri*,ri*) 



(63) 



k k' 
|k-k'| 



with Gvsikfirj,,) defined from G^siTi v) '''■^ in equation (43). The cancellation argument manifests itself here in 
that Gys{kiir]t,) is approximately the Fourier transform of a wide bell-shaped function Gys{r, t]o). If fcry* 3> 1, 
this implies that only for <^ 1 does Gys{k^r]i,) have strong contributions. Therefore we can set 7 ± k for 
the other slowly varying functions of /x. We are thus calculating the perpendicular part A^^ alone as the 
Vishniac term. 

We find that ^ 

Qvsik) = ^^r{kv*flr,,v6hMk)P\k,r].), (64) 
where Irj,vS is defined as in (61) in the obvious manner and 



Ik,vs{k) 



;i (1 - cos2g)(l ~ 2ycosef P[k{l + v^- 2y cosg)V^ 7?,] P[kyM 
dy / d(cos6') — — - — — — — -. (65) 



(l+y2_2t/cos6')2 



P{k,V*) 



P[k,ri^ 



For full ionization Ij^^ys — 6.38 when fig = 1 and ry* <C r]o- Again Efstathiou [6] uses this value for all cases. 
In the next section, we shall see what corrections are necessary for the general case. 

Care must be taken when evaluating the mode coupling [equation (65)]. The integrand has a divergence 
on the boundary at t/ = l,cos(? =1. As described in §IVA, a resonance occurs when the density field 
d^^^k') has a wavelength of A; ~ A;' (which implies y — 1) and the velocity field is coherent: |k — k'| <^ k 
{i.e., cos6 ~ 1). Furthermore, since k ± 7 to avoid cancellation, (k — k') || 7 and yields a strong Doppler 
contribution. The integral itself does not diverge however, since for reasonable power spectra, the power on 
the largest scales tends toward zero. Nonetheless, it complicates the integral, since we can have significant 
coupling between the largest and smallest scales. Thus even though for pure power laws Ik,vS is independent 
of k, we expect that realistic power spectra will deviate from this prediction. This is especially true for 
steeply red power spectra, where the large-scale velocity field is determined far above the scale of interest 
{i.e., k) where the power spectrum is quite different. Figure 1 displays the magnitude of this effect as a 
function of k for various choices of the small-scale power law n and an arbitrary large-scale cutofF in power at 
kmin = 0.001 Mpc~^. For reference the standard CDM model predicts n = —3, whereas BDM models prefer 
n = —1/2 but may plausibly lie in the range > n > — 1. This figure should only be taken heuristically 
since we have not included the proper large-scale power of a realistic power spectrum. 
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E. THE QUADRATIC DOPPLER EFFECT 

The quadratic effect escapes cancellation since it is not directly proportional to 7 • k = /i, but velocity 
perturbations are smaller than density perturbations at small scales, since in linear theory v oc 6/kr], if 
fio = 1 or {flozY^^ 3> 1. The effect is therefore relatively minimal at scales far below that of the horizon as 
long as the mode coupling is not significantly different from the Vishniac term. 

The source term for the quadratic Dopplcr effect is 

5,,(x,7, 7?) = -v(i) • v(i) + 7[7 • v(i)]2. (66) 
In Fourier space, a product becomes a convolution so that the solution can be written as 

A(2)(k,7,r?) = A,,(fcr?,)G,,(fcw*)exp(-ifcMr?)^F..(7,k,k')5(i)(k',r?,)^('Hk- k',r?0, (67) 

k' 

where 

A 



Gvv{r, v) 
F„„(7,k,k') 



2 

g{r],r'n*)r]^, (68) 



1 
2 

k' • (k - k') + 7(7 • k')7 • (k - k') fc^ 



Again G{k4ir]^,) is approximately the Fourier transform of a wide bell shaped function G{r,r]o), and thus we 
can set 7 _L k in the other terms. Notice that with this approximation 

_ fcfc^cos6> + k'^ - 7k'\m^ecos'^4) k^ 
~ fc2 + fc'2 - 2kk'cose A^' 

so the cross terms between this term and the Vishniac term, Fyg oc coscj), vanish under integration over the 
azimuthal angle (j). Furthermore, the quadratic Doppler effect produces spectral distortions in addition to 
anisotropics so that in principle they can be isolated from the linear Dopplcr effect. In practice however, we 
will find that the distortions from the quadratic Doppler effect are too small to be presently observable. 
The solution is 

1 y2 

Qvvik) = -^^ikr]:,)Ir,,vvIk,vvik)P'^{k,r]*), (69) 
TT-^ r]° 

where is defined as in (61) and 

'■+^ {y - cos0f - 7(1 - cos^ 9){y - cos9)y + (147/8) (1 - cos^Ofy'^ 



Ik,vv = dy d{cos9)- n , 2 o i 

Jo J-i (1 + 2/2 - 2j/cos6 

P[k{l + j/2 _ 2ycosey/^,r],] P[ky, r/,] 



(70) 



P[k,r]*] P[A;,r?* 



We have inserted the K = —2 factor from (32) under the assumption that we are observing the Rayleigh- 

Jeans temperature distortions from this effect. For full ionization, Qq = 1, and 77* ^ 779, we find Ir/.w — 3/2. 
In the general case, this value will change as does and Ir,^vS- This is because reflects the growth of 
perturbations as compared with their value at 77* , weighted by the probability of scattering at such an epoch. 
Perturbations grow more slowly in an open universe. If the integral is peaked above 77,, as is Lq,ys and Iri,vv, 
the integral will decrease; if it is peaked below 77*, it will increase, as with There is also a small effect 
due to cutting off the integral at the present time 770- Figure 2 shows the values of 7^ as a function of Os 
for Qo = and for flo = 1. 

Note that equation (70) is very similar to the Vishniac Ik,vS integral. Therefore the k dependence of 
Q{k) will resemble the Vishniac power spectrum save for the extra l/{kr]*)^ term suppression at small scales. 
Thus, unless 

^ > likrj^r^, (71) 
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wc expect the quadratic Dopplcr effect to be smaller than the Vishniac effc!ct at small angular scales. We 
calculate this ratio explicitly for a range of power laws, again with an arbitrary large-scale cut off in power 
at kmin = 0.001 Mpc~^ (see Fig. 3). In no case is the ratio sufficiently large to satisfy equation (71) for 
fcry* 3> 1. At scales nearer the horizon at last scattering, the linear first order Dopplcr effect is not cancelled, 
since there are only 0(1) scattering blobs of this size across the last scattering surface. Therefore, the 
quadratic Doppler effect plays no significant role at these scales either. 

IV. CALCULATIONS FOR SPECIFIC COSMOLOGICAL MODELS 



A. COLD DARK MATTER (CDM) SCENARIO 

Recent comparisons of anisotropy limits from the South Pole '91 experiment [3] with CDM normalized 
to COBE, have suggested that some CDM models, especially those with high Hb, may be ruled out, e.g., [25, 
26] . One possible resolution is to have early reionization to at least partially erase the primordial anisotropics 
on degree scales. If this were the case, then secondary anisotropics would tend to be generated on smaller 
scales, which we now consider in detail. 

The CDM power spectrum is given by [27] 

Pik, = Z — T2A7' (72) 

{l + [afc + (6A;)3/2 + (cA;)2]-}^/" 

where a = 6.4(r/i)-iMpc, b = 3.0{Th)-^Mpc, c = 1.7{Th)'^Mpc and u = 1.13. For standard CDM, F = 0.5. 
The normalization constant A is given by COBE to be [27] 

^ = 5.2xl0^f%Vi^^:^. (73) 

Figure 4 shows the radiation power spectrum Q{k) for the various effects in a universe that never recombines: 
Xe{z) = 1. For comparison, we have plotted the primordial fluctuations from standard recombination [28] in 
Fig. 4. As expected, first order and quadratic Doppler terms are suppressed at small scales. Furthermore, 
the quadratic effect never dominates on any scale. The result for the correlation function is plotted in Fig. 5. 
Here the correlations are given for infinite beam resolution a = 0. Note that the coherence scale is somewhat 
under an arcminute. 

To obtain observable quantities, the formulae of Wilson & Silk [29] can be employed, or we can note that 
for CDM on arcminute scales and greater, a double beam experiment essentially measures (7(0,0"), whereas 
a triple beam experiment measures 3/2 C(0, a) for the Vishniac effect, due to the small coherence scale. In 
Fig. 6, we show the predictions of C(0, a) for CDM with several values of Qb- 

Recently, Tcgmark & Silk [30] have shown that reionization may plausibly occur as early as Zi ~ 50. 
Furthermore, they show that reionization is sudden in realistic models, so that Xe{z) is essentially a step 
function at Zi. For these models, the optical depth will never reach unity until the standard epoch of 
recombination, and hence some primordial anisotropies remain. We will approximate this by employing the 
cumulative visibility function, 

rvo 

gciVi) = / 9{riQ,v')dri' (74) 



which tells us the fraction of the CMB photons which have been rescattercd since Zi and have thus suffered 
suppression of primordial fluctuations. In terms of the correlation function, only a fraction fprimiVi) = 
1 — gciVi) of the original signal remains after rescattering. 

Of course, the Vishniac effect decreases if the optical depth never reaches unity on the new last scattering 
surface. In this case, the epoch of new last scattering 77* is not the conformal time when optical depth reaches 
unity. We can however characterize this epoch by the max;imum of the visibility function. Employing this 
definition, we can express the fraction of the full effect for these models as fvs{i]i) = Qvs/Ql's'*^~^ ■ Figure 7 
plots the weighting fraction /(r/i) for the Vishniac effect and the primordial fluctuations. This fraction should 
be multiplied by the values in Fig. 4 to obtain the radiation power spectrum in cases when optical depth never 
reaches unity on the new surface of last scattering. Notice that we still obtain a fairly large Vishniac effect 
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even when only a small fraction of photons have scattered and the primordial fluctuations arc nearly intac;t. 
For example, if fls = 0.1 and Zi ~ 40 then fprim — fvS — 0.9, i.e., 90% of the primordial fluctuations remain 
and 90% of the maximum possible fluctuations from the Vishniac effect are generated. This is because the 
Vishniac effect becomes stronger at later times due to the growth of the velocities. Therefore even though 
only a small fraction of photons are rescattered a relatively large anisotropy is imprinted on them. So late 
reionization will not yield smaller CMB fluctuations on the arcminute level. 

We have also tested a phenomenologically successful model inspired by CDM that sets F — 0.2. This 
has the effect of suppressing the small-scale power with respect to the large leading to an extremely small 
Vishniac effect (see Table 1). 

B. BARYONIC DARK MATTER (BDM) SCENARIO 

Baryonic dark matter models with isocurvature fluctuations contain relatively large amounts of small- 
scale power (for the power spectrum in these models, which are also referred to as FIB or FBI models, see [31, 
32]. For this reason, primordial fluctuations from the first order Doppler effect at standard recombination 
are overproduced on the scale of the comoving particle horizon {i.e., rj). The angle subtended by the horizon 
at redshift z is 

At standard recombination z ~ 1100 and Oh — ■ Reionization is therefore necessary to remove these degree 
scale fluctuations. Generally, anisotropies will be erased up to the scale of the horizon size at the epoch z* 
where optical depth reaches unity. Secondary fluctuations from the first order Doppler effect on the new last 
scattering surface will of course be generated, and will peak near the horizon size at z*. Since the horizon 
scale is then larger than at standard recombination, the degree scale constraints on CMB fluctuations can 
be avoided. Primary fluctuations on degree scales will be erased and secondary fluctuations will be smaller 
because they have lower amplitude and because degree scales are well below their peak. 

In Fig. 8, we plot the angle subtended by the horizon at T(r?*) = 1. Notice that if f2o = 0.1 — 0.2 as 
suggested by some dynamical measurements, the horizon at last scattering subtends ~ 15° today. Therefore 
at degree scales, there is significant cancellation as described above. However, for the case in which the total 
density includes matter that does not contribute to the free electron density in the intergalactic medium 
(IGM), i.e., if baryons are hidden in compact objects or the ionization fraction is less than unity, then the 
last scattering surface is further away. In low universes, the angle subtended by the horizon at last 
scattering is only marginally larger than that at standard recombination. These models are therefore in 
danger of producing secondary fluctuations on degree scales which are as large as the primordial ones that 
have been erased. 

The relatively large amount of small-scale power in these models also leads to large, and in fact observ- 
able, fluctuations at arcminute scales due to the Vishniac effect. Thus even models with a relatively late 
last scattering epoch can be constrained. We normalize the spectrum to the fractional mass fluctuation on 
a scale of 8/i~^Mpc, i.e., as = 1. One complication arises however. We have performed our calculations in 
linear theory and therefore can only realistically predict fluctuations on scales larger than the non-linearity 
scale. We define the non-linearity scale as the value kni for which 

^ P{k,r,.)k'dk=l, (76) 

corresponding to the scale on which Sp/p ~ 1. The maximum anisotropy comes from taking the linear 
power spectrum out to infinitely small scales. The minimum anisotropy we expect comes from cutting 
off fluctuations at scales smaller than the non-linearity scale, i.e., setting Q{k > kni) ~ 0. Thus the 
true prediction lies somewhere between this minimum and maximum, the uncertainty being caused by our 
ignorance of nonlinear effects, although we would expect the prediction using the cutoff at kni to be more 
realistic. Choosing f2o = 0.2 and h = 0.7 as an example, we plot the results of using this prescription for 
various choices of the power law index n in Fig. 9. Note that for n > —1.5 the fluctuations due to the 
Vishniac effect increase with decreasing angular size, and therefore all experiments will effectively measure 
C7(0,a). 
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Recent measurements from the VLA have yielded significant improvements in the constraints on arc- 
second fluctuations in the CMB. Fomalont et al. [33] report C(0,o-)i/2 y. iq5 < ^.9, 2.1, 2.3, 4.0, 5.8, 7.2 for 
a = 0.56', 0.40', 0.22', 0.15', 0.10' respectively. Unfortunately, due to the uncertainty about effects below the 
non-linearity scale, the arcsecond measurements cannot be used to constrain most of the EDM models, and 
the larger angle measurements do not yet yield an improvement over the Readhead et al. [34] measurements 
at a = 0.78'. On the other hand, a small improvement in the measurements at cr > 0.2' would rule out many 
favored models, e.g., VIq = 0.2, h = 0.8, n = 0, x,, = 1 (see Table 1). 

The most stringent constraints to date come from the Australian Telescope Compact Array (ATCA) 
which uses interferometry to make a sky map and thus measure C(0, cr). Bounds are placed at u = 0.87' on 
C(0,(t)"'^/^ amounting to < 0.9 x 10^ [7]. This correspondingly places limits on the BDM model parameters 
fio and n (see Fig. 10). We have included a small contribution from the first order term Qy under the 
prescription that as ATCA does not measure structure above the primary beam size Opt = 2.5', there should 
be a Gaussian cutoff in the power above these scales. The tighter constraints on very low fio models are 
actually from this first order Doppler contribution. For these models, the last scattering surface is more 
distant, and the cancellation of the first order Doppler effect becomes correspondingly less complete. In 
general, whereas 0{v5) secondary anisotropics decrease with z,, 0{v) anisotropics actually increase with z*. 

Our constraints differ from Efstathiou [6], even accounting for the new limits. In particular, his scaling 
relations cannot be extended to certain regions of parameter space. The correction for the angle to distance 
relation i?,, [equation (35)] tends to increase the Vishniac effect for large n and decrease the first order term 
for small n. The true conformal time integrals 1^ (Fig. 2) tend to decrease the Vishniac term and increase 
the first order term for low CIq. The net effect is an increase in the predictions for low Oq and high n and 
a decrease in the predictions for low flo and low n. It is also interesting to note that the smaller scale 
experiments do a better job of constraining high n as opposed to low n. This is useful since the arcminute 
measurements tend to constrain low n (see Table 1). 

Notice that the ilo = 0.2, h = 0.8, n = — l/2,Xe = 1 model often quoted in the literature, e.g., [32], is 
ruled out by a small amount. We can avoid such constraints by lowering the density of free electrons in the 
intergalactic medium (IGM). An ad hoc procedure would be to lower the ionization fraction to a smaller 
constant value e.g., x,, = 0.1 [32]. But it is difficult to see how the ionization fraction can be kept at a small 
but significant value since the recombination time tends to be short. In these models, the predictions of the 
Vishniac effect for the ATCA experiment are then below the observational limits (see Table 1). On the other 
hand, we obtain an increase in the first order term. Again, this is because the last scattering surface becomes 
distant and cancellation is incomplete. If we decrease the ionization fraction substantially more than this, 
although primordial anisotropies will be erased, new ones of roughly the same size will be generated. These 
models may overproduce fluctuations on degree scales (see Fig. 8). 

A more physically motivated way of avoiding the constraint would be to hide the baryons in black holes 
[22]. Obviously, this has an effect quite similar to lowering the ionization fraction. For simplicity let us assume 
that the matter power spectrum is unchanged from the f^s = fio case [recall that we have defined fis as 
the fractional baryon density in the intergalactic medium (IGM)]. Fig. 11 shows how the combined Vishniac 
and first order anisotropies are changed by hiding a fraction of the baryons so that Q.b < ^o- The increase 
in fiuctuations for ^b/^q <C 1 is from the first order term and, as described above, is quite significant. For 
models in which the last scattering surface would already have been distant {e.g., f^o = 0.1, h = 0.5) the 
increase in the first order term is strong enough so that lowering the number of baryons in the IGM actually 
increases the predicted fluctuations on arcminute scales. Even in the favored f7o = 0.2, h = 0.8 model, hiding 
the baryons can only decrease the predictions by a factor of ~ 2. It seems that it may be possible to rule 
out even these models given modest future improvements in observations. 
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V. CONCLUSIONS 



We have shown that anisotropies on arcminute angular scales can provide a powerful probe of the 
ionization history of the universe. The COBE DMR detection of large-scale angular anisotropies implies 
intermediate angular scale anisotropies, (when combined with CDM or BDM power spectra) that in the 
absence of reionization, exceed or are comparable to observational limits. Reionization is essential in BDM 
models, and may be essential in CDM models, e.g., [25, 26]. Moreover, efficient reionization must have 
occurred at 2: > 5. Reionization has sufficiently modest energetic requirements (at least for photoionization) 
that it plausibly occurred when only a small fraction of the universe was found in nonlinear structures [30] . 
This occurs as early as 2; ~ 1000 in BDM and z 2± 50 in unbiased CDM models. While reionization helps 
suppress, or at least reduce, primordial fluctuations on degree scales due to the first order Doppler effect at 
standard recombination, it inevitably regenerates secondary fluctuations on arcminute scales. 

We have considered all Compton effects to second order and have shown explicitly that only those of 
0{v) and 0{v5) play a significant role in regeneration, i.e., no other higher order effects are important, and 
there arc no hidden terms to cancel the Vishniac term. Furthermore, we have substantially improved the 
approximations involved in calculating these effects in an open universe. The OixP') effect however does set 
the minimal spectral (Compton-y) distortions in reionization scenarios. For the CDM scenario, all of these 
secondary effects predict fluctuations well below present observations. However, one should note that even 
in cases where the primordial fluctuations are not efficiently erased due to relatively late reionization, the 
Vishniac effect still generates fluctuations that are of the same order as primordial fluctuations at arcminute 
scales. For BDM models where early reionization is necessary, new measurements from ATCA already 
exclude a signiflcant region of (^lo,?^) parameter space. In fact, the standard model [32] of VIq — 0.2, h = 0.8 
with full ionization = 1 is ruled out by a small amount. Changing the ionization history can avoid the 
present constraints, but it is difficult to decrease the secondary fiuctuations by a significant amount. If 
one arranges the ionization history to produce a small Vishniac effect, the first order Doppler effect will be 
correspondingly increased. We therefore find for both CDM and BDM that changing the ionization history 
in what might be considered physically well-motivated prescriptions does not in fact greatly change the 
predictions for the secondary anisotropy. 
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APPENDIX A: COLLISIONLESS BOLTZMANN EQUATION 

Second order gravitational effects play an important role when the first order gravitational effect vanishes, 
e.g., if there are no metric fluctuations on the last scattering surface as in the case of late phase transition 
scenarios. These effects will primarily manifest themselves on large scales. Martinez-Gonzalez, et al. [35] 
have performed the calculation for weakly non- linear effects due to second order density perturbations. Here 
we complete the program by deriving all possible effects to second order in the metric perturbations. Note 
however that recently Jaffe et al. [36] have shown that scenarios in which there are no fluctuations on the last 
scattering surface must involve strongly non-linear effects to generate the observed gravitational potentials 
today. They find that such models must generate CMB fluctuations on the order of, or greater than, the 
standard primordial fluctuations. Therefore, even for this class of models, weakly nonlinear effects such as 
those derived in here may only play the role of a small correction to the full gravitational effect. 

We choose to work in the synchronous gauge, following closely the notation of Peebles and Yu [37] . In 
this gauge, perturbations to the metric are expressed as 

300 = 1, 5oi = 0, gij = -a{tf[Sij - hij{x,t)], (A-1) 



where greek indices run from 0-3 and latin indices run from 1-3. Summation over repeated indices is assumed 
throughout even when all indices are lowered. The momentum of the photon will be denoted by and for 
convenience we express, 

Pi = -poa{t)e'yi, (A-2) 
where 7^ are the direction cosines. In order to satisfy p^p^ = 0, 

= 1 - hij-yaj, (A-3) 

to first order. 

The equations of motion for the photons are [37] 

-TT = i;9agi^^ . (A-4) 

dt 2 Po 

The general expression for the coUisionless Boltzmann equation is 

dt dx' dt d'ji dt dpo dt ' ^ ^ 

We integrate this equation over energy to obtain the equation for A, the brightness. Each term in the 
integral of (A-5) is now discussed in turn. 

1. Redshift Term: 



t J opo dt a dt at at (A-6) 



dhij 



dt ' 

where £ is the average energy density of the photons. The first term represents the universal redshift due to 
the expansion. Note that the brightness fluctuation is independent of the redshift factor since all photons 
redshift in the same manner. Thus, the term from the explicit time dependence of the spectrum will cancel 
the first term (see 4.). The other terms represent gravitational redshifts due to perturbations of the metric. 
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2. Anisotropy Term: 

Anisotropy is generated as a first order gravitational effect due to gravitational redshifts from overdense 
regions, i.e., the Sachs-Wolfe effect [38]. Therefore df /dji has a contribution to first order in the metric 
perturbations. Gravitational lensing, d-yi/dt, is also a first order effect. Thus the whole term is second order: 

47r /■ 3 , df d-y, 5A(i) 1 r dh.k 1 o , 1 , , „^ 

Y I ^0*° d^,^ = YW^'^^^' + -d,hm] . (A-7) 



3. Inhomogeneity Term : 



47r /• , df dx' r aA(i) SA^^) 1 l aA(i) 



£ J dx^ dt a 1 dx^ dx^ ( a dx^ 



(A- 



This term represents the effects of a spatially dependent perturbation in the photon energy density. The 
first term has no dependence on the metric fiuctuations and is the fiat space approximation used in equation 
(22). The second term represents corrections due to the mass shell constraint, equations (A-2, A-3). 

4. Explicit Time Dependence Term : 



_ y ,3,^^ ^ ^ ^ ^(r) ^ ^(.)) + ^ + (A-9) 

The first term represents the uniform redshift of the spectrum. It cancels with the energy redshift leaving 
no effect on temperature perturbations. 

Complete First Order Expression: 

Summing the above, we arrive at 

to first order. This term represents the gravitational redshift of the photon due to the metric fiuctuations 
and gives the conventional Sachs- Wolfe effect [39]. 

Complete Second Order Expression: 

In second order. 



5A(2) 7,aA(2) dhf/ r .dhij 



aA(i) 1 r dhik 1 o , 1 If 9A(i 



d^i 2 t Si •"'^ ■ a " '"J a t dx^ 

(A-11) 

This expression gives the second order Sachs- Wolfe effect and gravitational lensing effects. 
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APPENDIX B: COMPTON COLLISION TERM 

Here we present the general expansion of equation (3) to second order, and give the exphcit sources in 
the colhsion term of equation (6). First we must express the matrix element in the frame of the radiation 
background: 



|M|2 = {^^^0^21 (1 + cos^/3) - 2cos/J(l - cos/3) 



q • p ^ q • p 



mp mp 



+ (l-cos/3)2^ 



+ ^cos/3(l - cos/3) + (1 - cos/3)(l - 3cos/3) 



q p ^ q p 

mp mp' 



2cos/3(l - cos/3) 



(q-p)(q-p') 



m^pp' 



h.o., 



(B-1) 



where higher order (h.o.) represents third and higher order in v and p/m. The following identities are very 

useful for the calculation. Expansion to second order in energy transfer can be handled in a quite compact 
way by denoting it as an expansion in 6p = q — q' of the delta function for energy conservation: 



1 



Sip + q-p'-q')=S{p-p') + — {2{p-p')■^+{p-p'f} 



dp' 



-6{p-p') 



+ g;^{2(p-p')-q+(p-p')'}' 



dp 



5{p-p') 



(B-2) 



+ h.o. 



which is of course defined and justified by integration by parts. Integrals over the electron distribution 
function are trivial, 

d^q , , 



(27r)3 



I 



(B-3) 



Now we substitute all this into the general collision equation (3) and obtain the explicit expressions for the 
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terms of equation (6), 



Co =d{p - p') (1 + cos2/3) i^i(x, p, p') 



Cv 



a 



p/v 



a 



Te/m 



(1 + cos"/3) V • (p - p') 

V • p V • p' 



1 


52 


2 













-(5(j3-p')2cos/?(l - cos/?) 



l + cos2/3) [v(p-p')]'Fi(x,p,p') 



i^i(x,p,p') 



-5{p-p') 



2cos/3(l - cos/3) 



V • p V • p 

— i: + 

P 



p' 



+ S{p -p')< - (1 - 2COS/3 + 3cos2/3)t;2 + 2cos/3(l - cos/3) 



v(p-p')Fi(x,p,p') 
(vp)(vp') 



pp 



+ (1 - cos/3)(l - 3cos/3) 



V • p V • p 



p p' 



>i^i(x,p,p') 



(1 + cos^/j) 



(P-P') 



2 m 



-^2(x,p,p') 



(l + cos^/3)iP^ 



/^2 



2cos/3(l-cos2/3)(j9-p') 



+ S{p - p')[4cos^p - 9cos2/3 - 1] >:^Fi{^, P, p'). 



(B-4) 
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and the higher order terms, 

+ 
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(1 + cos^/?) (p - pf ^-(P^ P'^ f,(x,p,pO 



2cos/3(l -cos^/3)(p-p') 



2m 

,^v.(p-p') 



m 



^3(x,p,p') 



+ cos/?(l - cos/3)(p - p')' ( ^ + ^ ) F2(x, p, p') I 

\ mp mp' J I 

+ S{p - p') I ~ [- (1 + cos^/?) Fi (x, p, p') + 2(1 - 2COS/3 + 3cos2/3)F3(x, p, p')] 



2cos/?(l - cos/?) 



P fl\ v-p' 

cosp 

p' J m 



I /ml 



cosp 



m 



i^3(x,p,p') 



[p/my — 
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ifi(x,p,p') + F3(x,p,p') 



2cos/3(l - cos2/?)(p -p')-^^2^ 



i^3(x,p,p') 



+ -cos/3)2 (^)'fi(x,p,p') + (l + cos2/3) IH^-^F^Ix, p, p') 
- 2cos/3(l - cos^)^ F3(x,p,p') 



- [1 - 2COS/3 + 3cos^/3] 



(p-p'r 



F3(x,p,p') 



(B-5) 



where we have written 



for compactness. 



Fi(x,p,p') =/(x,p')-/(x,p) 

F2(x, p, p') =/(x, p) + 2/(x, p)/(x, p') + /(x, p') 

F3(x,p,p') =/(x,p')[l + /(x,p)] 



(B-6) 
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FIGURE CAPTIONS 

FIG-1. Mode coupling integral for the Vishniac effect {vS) [see equation (65)] for various power laws [P{k) oc 
fc" for large k] with a large-scale cutoff at kmin = 0.001. Note that for steep power spectra the integral never 
approaches a constant. This represents a correction to the approximations in ref. [6]. 

FIG-2. The conformal time integrals [sec equation (61)] for the 0{v), 0{vS), and 0{vv) effects for various 
models. Note that for 1 the values depart significantly from the asymptotic = I, f^s <C 1 result. 

Again this significantly improves the approximations of ref. [6] . 

FIG-3. Ratio of mode coupling integrals, i.e., Ik of the quadratic Doppler effect (vv) to that of the Vishniac 
effect (vS). For no power law does the ratio become sufficiently large to make the quadratic Doppler term 
significant compared with the Vishniac term, at small scales. 

FIG-4. CDM radiation fluctuation power spectrum for the first order Doppler (v), Vishniac {v6), and 
quadratic Doppler (vv) effects. Here Q(fc) is the power per logarithmic interval, i.e., C{0) — J Q{k) dk/k. 
The shape parameter F [see equation (72)] is 0.5 for standard CDM. Note that the {vv) term is never 
significant: (v) dominates at large scales and {v5) dominates at small scales. 

FIG-5. CDM predictions for a double beam experiment with infinite resolution. The primordial fluctuations 
for standard recombination are normalized to COBE [2]. Note that at sub-arcminute scales the Vishniac 
(vS) term dominates over the first order Doppler term (v) and can be larger than the primordial fluctuations. 
Also note that the coherence scale of the Vishniac effect for CDM is somewhat under an arcminute. 

FIG-6. CDM predictions of the Vishniac effect for C(0, a) for various choices of Hb- Because the coherence 
scale is under an arcminute, experiments at arcminute scales and greater will essentially only measure C(0, a). 

FIG-7. Weighting fraction for CDM models with sudden reionization at Zi. The full values of Figure 4 should 
be weighted by / to determine fluctuations for these ionization histories. Notice that we still obtain a strong 
Vishniac effect even in cases where the primordial fluctuations are only partially erased. At Zi ~ 40, we have 
approximately 90% of both the Vishniac and primordial fluctuations. Late reionization therefore does not 
necessarily yield smaller fluctuations. 

FIG-8. Angle subtended by the horizon at last scattering (optical depth unity), i.e., rj^, [see equation (75)]. 
In models where fis = ^o, the horizon size is sufficiently large to escape degree scale constraints. However, 
models which hide matter in compact objects or have low ionization fraction will have horizon sizes close 
to the degree scale. We have plotted for illustrative purposes a model in which 90% of the baryons are in 
compact objects (fis = O.IOq) or alternatively ionized fraction Xg = 0.1. 

FIG-9. EDM predictions for C(0, a) from the Vishniac effect with Oq = 0.7, h = 0.7, Xe = 1.0. This model 
barely escapes constraints by ATCA. Notice that at arcseconds we are in the non-linear regime; unccat,aintics 
here greatly inhibit our ability to constrain models. We have plotted predicted fluctuations with (heavy) 
and without (light) the cutoff at kni- These should be taken as lower and upper bounds. 

FIG-10. The limit from the ATCA measurement of C(0, 0.87') < 0.9 x 10~^ places signiflcant constraints on 

BDM models with full ionization Xe = 1. The constraints are strong fimctions of the Hubble constant; here 
we show three representative values of h, the lowest and highest values generally considered and the specific 
value {h = 0.8) chosen by ref. [32]. The boxed region indicates the preferred values of ref. [32]. Low n and 
high rig are excluded. Note that other assumptions about Xe would alter these limits. 

FIG-1 1. Combined first order and Vishniac predictions for BDM models with a significant fraction of the 
baryons in compact objects which do not contribute to the baryon density in the intergalactic medium, Qb- 
The first order effect dominates for low CIb, since the last scattering surface is then at such high redshift. 
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TABLE 1 



Model 




h 


spectrum 




a = 0.1' 


= 0.5' 


= 1.0' 


BDM 


0.2 


0.8 


n = 


1.0 


2.61-11.50 


1.17-1.19 


0.48 


BDM 


0.1 


0.8 


n = -0.5 


1.0 


2.17-3.26 


0.80 


0.48 


BDM 


0.1 


0.5 


n = -0.5 


1.0 


1.77-2.36 


0.49 


0.25 


BDM 


0.2 


0.8 


n = 


0.1 


1.23-1.35 


0.13 


0.05 


BDM 


0.2 


0.8 


n = -0.5 


0.1 


0.83 


0.18 


0.09 


CDM 


0.1 


0.5 


r = 0.5 


1.0 


0.38 


0.25 


0.18 


CDM 


0.1 


0.5 


r = 0.2 


1.0 


0.04 


0.03 


0.03 



TABLE CAPTION 

TABLE-1. Predietions of C(0, cr) (xlO~^) for the Vishniac effect for various models which can be tested 
by future experiments at the VLA, ATCA, and OVRO. We have chosen three representative values of the 
beam smoothing a. Ranges given in the 6th and 7th column represent the uncertainty in predictions due to 
nonlinearity. The observed signal will also depend on the first order effect which in turn depends sensitively 
on the nature of the experiment and not on C'{0,a) alone. 
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